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In a new branch of quantum computing, information is encoded into coherent states, the primary 
carriers of optical communication. To exploit it, quantum bits of these coherent states are needed, 
but it is notoriously hard to make superpositions of such continuous- variable states. We have realized 
the complete engineering and characterization of a qubit of two optical continuous- variable states. 
Using squeezed vacuum as a resource and a special photon subtraction technique, we could with 
high precision prepare an arbitrary superposition of squeezed vacuum and a squeezed single photon. 
This could lead the way to demonstrations of coherent state quantum computing. 



Among the various physical implementation schemes 
of quantum information processing (QIP), optical QIP in 
traveling light fields is a significant contender [1 j. Used 
at the nodes of a quantum optical network, it would de- 
liver ultra-high capacity with minimum power in opti- 
cal communications [2] and highly secure communica- 
tions and distributed QIP [3H5] surpassing any classical 
counterpart. Unfortunately, photons do not readily in- 
teract with each other, making it quite a task to im- 
plement quantum gate elements. One currently feasible 
scheme is linear optical quantum computing (LOQC), 
which uses off-line resource states, linear optical process- 
ing, and photon- number resolving detection [6] [7]. There 
are two approaches to LOQC, the standard one being the 
single photon scheme, where single photons are used as 
the physical quantum bits (qubits) [6 . The other is re- 
ferred to as coherent state quantum computing (CSQC), 
where two phase-opposite coherent states are used for the 
qubits, i.e. |t) = \a), \\) = \-a) [SHE]. 

CSQC is not only effective for exponential speed-up 
of computations, but also for attaining the ultimate ca- 
pacity of an optical channel in current network infras- 
tructure where coherent states are the primary carriers. 
Because coherent states propagate intact, even through 
lossy channels, simple coherent state-encoding is found to 
be the optimal transmission strategy [T3] , At the same 
time, the optimal decoding should be fully quantum, and 
can be implemented by an extension of CSQC [14]. Prac- 
tical implementation of CSQC is still a big challenge, 
though - one requirement is the availability of arbitrary 
qubit states as resources. So far, two diagonal states of 
the qubit \a) ± |— a) - so-called Schrodinger cat states - 
have been generated in the laboratories [T5HT9] . 

We have implemented a setup that is suited for the 
generation of such arbitrary qubits. For this demonstra- 
tion, though, we perform the complete engineering of a 
different, but closely related kind of qubit, namely one 
with squeezed vacuum and squeezed single-photon states 
as the basis. The squeezed photon state is in fact very 
similar to one of the CSQC diagonal qubits, as was uti- 
lized in the previous demonstrations of those. To create 



the arbitrary superposition of the basis states, we use 
single-photon subtraction from a squeezed vacuum as- 
sisted by a coherent displacement operation [20] . Photon 
subtraction is a simple, but powerful technique that has 
been used for non-classical state generation [T5UT9] . en- 
tanglement increase [2TJ [22], and fundamental tests of 
quantum mechanics [23] (see [24] for an overview). A re- 
lated experiment for engineering of superpositions of the 
0-, 1-, and 2-photon states was recently reported [25] . 
In our scheme, in contrast, the state elements of the su- 
perposition lie in the infinite dimensional Hilbert space, 
including many-photon number states. We should also 
note that spectacular progress in generation of complex, 
high-photon number states has been made in microwave 
resonator fields [26j |27] , but these states are trapped and 
therefore of limited usefulness. 

A simplified schematic is shown in Fig. [TJl. The 
squeezed vacuum state Sa (r) \0a) is prepared in mode 
A, where Sa (/*) represents the squeezing operation with 
strength r. A small fraction R of it is tapped off via a 
beam splitter as a trigger beam in mode B, subjected to 
the displacement operation Db(P), and detected on an 
avalanche photodiode (APD). Without the displacement 
operation, the output state conditioned on a click at the 
APD would be 

(1b\ VabSa(t) \0 A b) « -VRsmhrS A (r) \1 A ) , (1) 

where Vab is the beam splitting operator. That is, a pho- 
ton has been subtracted from the squeezed vacuum state 
in mode A, which is equivalent to squeezing a single pho- 
ton. With the inclusion of the displacement operation, 
this changes to 

(l B \DB(P)VABS A (r)\0AB) « 

Af(l3S A (r) \0 A ) - VRsmhrS A (r) \1 A )), (2) 

with normalization factor Af. This superposition origi- 
nates from the two indistinguishable processes; an APD 
click comes either from the displacement or from the 
squeezing, corresponding to the output states S(r) |0) or 
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FIG. 1. Experiment outline, a) Conceptual schematic, 
b) Experimental setup. OPO, optical parametric oscilla- 
tor; R, beam splitter reflectivity; LO, local oscillator; APD, 
avalanche photo diode. 



S(r) |1), respectively. A displacement photon is uncorre- 
cted with the output mode, hence leaving the squeezed 
vacuum state intact. The superposition weight and phase 
can be completely controlled by the displacement opera- 
tion. Each of the two states is composed of several Fock 
state elements - only even photon numbers for S(r) |0) 
and odd numbers for S(r) |1). They are orthogonal to 
each other, so together they constitute a qubit basis and 
the general output state ([2| can therefore be represented 
on a Bloch sphere as 

\p(6, ^ = cos 6 -S(r) |0> + e* sin °-S{r) |1) , (3) 

with (p = 7r — arg f3 and, after normalization, 
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where ndi sp (^sq) are the number of photons in mode B 
originating from the displacement beam (squeezing). 

Our experimental setup is shown in detail in Fig. ^p. 
The input squeezed vacuum states are in a cw beam gen- 
erated by an optical parametric oscillator (OPO) at a 
center wavelength of 860nm with a bandwidth (o/tt ~ 
9 MHz (FWHM). This cw beam is intuitively a continu- 
ous sequence of squeezed light packets, each of which is 
in a temporal form ij)(t) = ^/Coe~^ ^ (the squeezing level 
within these packets is -2.7 dB in our experiment). 

After tapping off R = 5% for the trigger, the main 
part of the light (the output signal) is measured on a 



homodyne detector for state analysis. The trigger beam 
is spectrally filtered by two subsequent Fabry-Perot res- 
onators before it is displaced and directed onto a Si APD. 
The phase space displacement is implemented by overlap- 
ping the beam with a weak coherent state on an imbal- 
anced beam splitter [28]. For experimental convenience, 
we use a combination of half-wave plates (HWP) and 
a polarizing beam splitter (PBS) which allows for inde- 
pendent tuning of the splitting ratios of the two beams. 
For the trigger beam, 90% is reflected towards the APD, 
while only ~0.1% of the displacement beam is transmit- 
ted in that direction. The other output mode of the PBS 
is monitored on a standard linear photodiode for the pur- 
pose of locking the displacement phase. The phase mon- 
itoring is done with the help of chopped probe light (10 
kHz, 20% duty cycle) in both the OPO output and the 
displacement beam modes, and the locking is controlled 
by FPGA modules. Apart from this essential displace- 
ment part, the setup is mostly identical to that in Refs. 
[17j[T9], where more details are provided. 

The time resolution of the trigger signal is on a scale of 
sub ns, which is almost instantaneous compared with the 
time scale of the squeezed packets (tt/(o ~ 100 ns). The 
trigger signal, say at time t±, specifies a temporal mode 
ip(t — ti) of the output state of interest. In the homodyne 
channel, a continuous photo-current is sampled around 
the trigger time t±, and subsequently integrated over this 
packet mode function to yield the observed quadrature 
variable of the field. The specific quadrature to be mea- 
sured is determined by the phase of the local oscillator 
(LO) of the homodyne detector. To obtain full informa- 
tion about the output quantum state, we must carry out 
many homodyne measurements at a range of different LO 
phases, that is, a tomographic state reconstruction. For 
each state, 360,000 quadrature points were observed, dis- 
tributed on 12 different LO phases. The reconstruction 
was then done by the maximum likelihood method [29] 
without any correction of measurement losses. 

A selection of our generated states are presented in 
Figs. [2] and [3J Each state is represented by its Wigner 
function as a top-down contour plot and a 3D view, as 
well as by a Bloch sphere map of the fidelity between 
the state and the ideal squeezed qubit ([3| for all combi- 
nations of and (j). To demonstrate the performance of 
the state engineering, we show in Fig. [2] the control of 
the superposition weight while keeping the phase con- 
stant at a) (f = 0° and b) ip = —90°, and conversely, in 
Fig. [3] we show the control of the complex phase cp for 
fixed weights of the superposition, a) = 135° and b) 
= 60°. In Fig. [2] we see that by increasing the amount 
of displacement in the trigger channel, |/?|, we can move 
from the south pole (squeezed photon) to the north pole 
(squeezed vacuum) of the Bloch sphere along a fixed lon- 
gitude. While doing that, the negative dip of the Wigner 
function moves away from the center in a direction de- 
termined by the displacement phase. At the same time, 
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FIG. 2. Control of superposition weight. A variety of experimentally generated superposition states where the superpo- 
sition weight has been swept with the phase <j) fixed at a) 0°, and b) —90°. The left and center panels show the Wigner 
function as a contour plot and surface plot, respectively. The axes are the x- and p-quadratures. The right panels show a 
flattened Bloch sphere with an overlay signifying the fidelity between an ideal qubit (eq. [3| at the given cp) point and the 
measured state. The central longitude corresponds to (p = 0°, and = 0° at the north pole. The small circles serve to point 
out which state we were aiming at. The ideal qubit is taken to have a squeezing parameter r = 0.38. 



the dip becomes shallower and finally disappears as the 
state approaches the north pole. In Fig. [3j on the other 
hand, when sweeping the displacement phase, arg/3, the 
negative dip circles around the center while the state on 
the Bloch sphere turns around at a fixed latitude. 

From the fidelity maps, we can see that there is a 
clearly defined qubit state of maximum fidelity in the 
center of the orange parts of the maps. These maxi- 
mum points are all quite close to the target states that 
we aimed for in each qubit realization - these targets 
are marked by small circles on the fidelity maps. This 
illustrates the precision of our state control. It is also 
clear that the obtained fidelities are not as high around 
the south pole as they are in the north. That is be- 
cause highly non-classical states with negative Wigner 
functions, such as the squeezed single photon, are much 
more fragile and susceptible to losses than, for example, 



squeezed vacuum. 

Equation Q anticipates a direct correspondence be- 
tween the APD click rates and the resulting superposi- 
tion weight. The photon number from squeezing (dis- 
placement) is directly proportional to the count rate R sq 
(^disp) observed when blocking the displacement (trig- 
ger) beam before the overlapping PBS, so the relation 
would read = 2 tan _1 ( J R c iisp/^sq) _1 ^ 2 - This relation is 
shown as the green curve in Fig. [4j where also the exper- 
imentally obtained (of the ideal qubit with maximum 
fidelity) as a function of the click rate ratios i?disp/^sq are 
plotted. The data points are shown for both the (p = 0° 
and the ip = —90° series of state generation. All the 
points are lying below the theoretical curve. That is be- 
cause the model in eqs. ([2j|4]) is an idealized picture from 
which there are several deviations in the actual experi- 
ment. In particular, the unavoidable losses in state gen- 



FIG. 3. Control of superposition phase. Generated states with the superposition phase (f) swept while keeping the 
superposition weight constant at a) 135°, and b) 60°. 



eration and measurement serve to mix in vacuum which 
effectively acts to decrease the obtained 0. To more ac- 
curately describe the experimental outcomes, we have 
developed a relatively simple model (see the appendix) 
which takes into account a more detailed description of 
the OPO and the trigger filtering, as well as all the ineffi- 
ciencies of the setup. This model, with no free parameters 
(except for the squeezing parameter r that is semi- fixed), 
is also plotted in Fig. [4] and is seen to simulate the mea- 
sured outcomes very well. Apart from the superposition 
weights, the figure also shows the fidelities between the 
measured states and the target states - that is, the fi- 
delity values at the target marks in the Bloch sphere 
maps of Fig. [2] The state preparation works somewhat 
better for displacement along the anti-squeezing direction 
((f) = 0) than along the squeezing direction. This can be 
ascribed to the fact that losses have a larger influence 
on the squeezed than on the anti-squeezed quadrature. 
The fidelities are also well described by the theoretical 



model, although with some smaller discrepancies in the 
weak-displacement regime. 

As these results show, we were able to realize with high 
precision and relatively high fidelity the complete engi- 
neering of a qubit of continuous- variable states. The sim- 
ple superposition preparation technique demonstrated 
here can be straight-forwardly applied to generation of 
coherent state qubits. If the input state instead of 
squeezed vacuum were one of the 'cat' states \a) ± |— a), 
then the single photon subtraction would change this in- 
put into the opposite state \a) =p |— a). Thus, combining 
this technique with the already existing cat state [19] , we 
should be able to generate an arbitrary superposition of 
these two cat states. This is equivalent to arbitrary su- 
perpositions of the coherent states \a) and \—a), i.e. the 
coherent state qubits that are cornerstones of CSQC. As 
a side-note, the currently measured basis states, squeezed 
vacuum and squeezed photon, already have a quite good 
resemblance with the '+'-cat and '-'-cat states with fideli- 
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FIG. 4. Influence of displacement strength. Experi- 
mentally obtained superposition weights (filled points + solid 
curves, left axis) and fidelities with the intended target states 
(unfilled points + dashed curves, right axis) versus the APD 
count rate of displacement photons relative to squeezing pho- 
tons for two series of measurements with different superposi- 
tion phase factor. The curves are obtained from a theoretical 
model with no free parameters, taking into account the var- 
ious imperfections. The green curve is the relation for the 
idealized, loss-less case. Note that the extreme data points 
correspond to count rate ratios of and infinity. 



ties of 81% and 68%, respectively, for a cat-amplitude 
a = 1.0. It is also possible to generate larger ampli- 
tude cat states by cascading the photon-subtraction [30] . 
The results reported here can therefore lead the way to 
prototypes of CSQC quantum gates, likely to become im- 
portant ingredients for attaining the ultimate capacity in 
future quantum optical networks. 
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Appendix: Theoretical model 

In this appendix, we outline the theoretical modeling 
of the quantum states generated in our experiment - a 
model that fits very well with the observations, as seen 
in Fig. 4 of the main paper. 

We start by recounting the experimental setup and lay- 
ing out all the relevant parameters. Next, we describe 
a general Wigner function formalism for calculating the 
conditional states and relate it to our setting, including 
imperfections of the displaced photon subtraction. Fi- 
nally, we specify the actual input Gaussian state, deter- 
mined by the OPO output correlations, the trigger chan- 
nel filtering, and our choice of temporal mode functions. 
To round off, we also include the Wigner function for the 
ideal squeezed vacuum/squeezed photon qubit - the class 
of states that we match fidelities against. 

Other relevant and more general theoretical descrip- 
tions can be found in Refs. [IH3]. 



EXPERIMENTAL SETTING 

Initially, a cw squeezed vacuum is generated by a sub- 
threshold OPO into mode A. The OPO HWHM band- 
width is 7 and the pump level is e < 7. We align the 
phase space to have squeezing in the p-quadrature. The 
normal ordered temporal output correlations are then [4 



(:A£ A ,mi(t)A£ A ,ini(0:) 



^ c -fr-e)\t-t'\ 

7-e 



(:Ap A ,ini(t)ApA,ini(t ):) 



'V\ ^ r -h+e)\t-t'\ 



7 + e 



(5) 
(6) 



The squeezed vacuum is split on a beam splitter with 
transmission T £ , being mixed with vacuum from mode B. 
The transmitted beam is the output signal to be recorded 
by the homodyne detector whose overall efficiency we de- 
note by tja ~ this efficiency parameter includes internal 
OPO losses, propagation losses etc. The reflection is sent 
towards an APD, being displaced on the way (after fre- 
quency filtering) by mixing it with a coherent beam. The 
details of the displacement process are not important: 
What matters is the rate of APD clicks originating from 
the squeezed light, i? sq , or from the displacement beam, 
Rdisp ~ a click coming from the squeezed light heralds a 
photon subtraction in mode A, while a click coming from 
the displacement beam (uncorrelated with mode A) her- 
alds no action. 

Other important parameters are the spatial mode 
matching between trigger and displacement beam, x, the 
dark count rate, R dc , the overall efficiency of mode B 
(trigger), r] Bl the trigger filter bandwidth [5], ft, and the 
phase angle of displacement, </>di S p- The spatial mode 
matching x must be close to 1. If not, the two possible 
origins of the detected photon are not indistinguishable, 
leading to a mixed state instead of a superposition. 



The parameters used in this model together with their 
typical experimental values are listed in Table [IJ Apart 
from these, there are also the parameters that define the 
signal mode function as used in the extraction of the 
homodyne data: 7/, e/, and Kf. These can be chosen 
freely, but we usually use the experimentally-equivalent 
values, that is, 7/ = 7, €f = e, Kf = k. 

GENERAL OUTPUT STATE FORMALISM 

To describe the effect of the photon subtraction, we 
adopt the Wigner formalism of Ref. [6] . The two- mode 
state prior to the photon detection event is Gaussian, i.e. 
it has a Gaussian Wigner function fully determined by 
the covariance matrix T and displacement vector d: 



W G (x A ,p A ,x B ,p B ) = 



1 



yrVdetr 



-(x-d) T r- 1 (x-d) 



(7) 



with x = (x A ,p A ,x B ,p B ) T . 

To get the state conditional on an APD click in mode 
B, the Gaussian Wigner function should be multiplied 
by the Wigner function corresponding to the detection 
operation and integrated over mode B - this corresponds 
to ti B [p AB Il B ] in the density matrix formulation. For 
the non-photon number resolving APD, we use the stan- 
dard on/off operation II on = I — |0) (0|, which has the 
equivalent Wigner function 



W on (x,p) W [d (x : p) - W vac (x,p) 



1 1 

2tt 7T 



-x —p 



(8) 

The output state conditioned on a displaced 1-photon 
subtraction is therefore 

W d -i ps (x A ,PA) 

M out 27ij I W G (x A ,p A ,x B ,p B )W on (x B ,p B ) dx B dp B , 

(9) 

with normalization to be determined by integrating the 
output state. 

Including dark counts and mode matching 

If an APD click is a dark count or if it comes from the 
part of the displacement beam that is not mode-matched 
to the trigger beam, there is no action on the signal beam 
and the output state is simply a squeezed vacuum state 

W sq (x A ,p A ) = 

2tt W G (x A ,p A ,x B ,p B )W id (x B ,p B )dx B dp B 
J J —00 

•00 

W G (x A ,p A ,x B ,p B )dx B dp B . (10) 
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OPO bandwidth 


2ttx4.5 MHz 


e 


pump level 


variable - in this paper, O.37 


K 


filter bandwidth 


- 2ttx25 MHz 


T t 


tapping BS transmission 


0.95 




overall efficiency, signal 


0.82 (0.96 x 0.91 x 0.98 2 x 0.99 x 0.99) 


VB 


overall efficiency, trigger 


- 0.1 


Rsq 


click rate, squeezed photons 


variable - in this paper, 3600 c/s 


Rdc 


click rate, displacement photons 


-30 c/s 


-Rdisp 


click rate, dark counts 


variable - typically 50-50,000 c/s 


0disp 


displacement angle 


variable 


X 


displacement mode matching 


0.97 



TABLE I. Model parameters and typical experimental values. The factors of the signal efficiency are: OPO escape, propagation, 
LO visibility, quantum efficiency, dark noise. 



If, on the other hand, the click originates from the trigger 
beam, but from the part which is not mode-matched to 
the displacement beam, the result will be a normal un- 
displaced photon subtraction: 

Wi ps (x A ,p A ) 

Mout^l I WG, n &(xA,PA,XB,PB)W n{xB,PB)dx B dpB- 
J J —00 

(ii) 



with the un-displaced Gaussian state 

W gm (x a ,Pa,x b ,Pb) = ~ 1— - e~ xTr ~ lx . (12) 

7r 2 vdet T 

Let the total count rate be R = R sq + i^disp + ^dc- 
The final output state is then (with (xa,Pa) temporarily 
removed) 



Wont — 



disp 



R 



xWd-ip S + (1 - x) 



R 



disp 



^sq + Rdh 



-w sc 



R s 



Rs 



R 



sq 



Rd] 



disp 



R 



w d _ 



R 



lps 



+ (i-x) 1 fWV+((i-x) 



v ^di 



sp 



R 



Rdisp 

Rd 
R 



~R W ^ 



(13) 
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Explicit Wigner function 

Any two-mode Gaussian state can - via local opera- 
tions - be put on the following generic form: 



(14) 



By having aligned the squeezing quadrature of the initial 
state along one of the phase space axes, we have already 
obtained this form (see next section), where the x and p 
variables are completely decoupled. We can set r = s = 
as displacement only occurs in mode B. 

When inserting Y and d in ^ and carrying out the 
integrations, we get the following normalized expressions 





/ e ^ 
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for the three state components of (13): 
1 



W sq (x Aj p A ) 

Wi ps (x A ,p A ) 
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1-2 12 
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l-w Tr^/ab 
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1 1 _ I r 2 _ 12 
-g a X A bPA 
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1 - w d ttVVF 



-e a 



jjixA-rdf-^ripA-SdY 



with 



a = a 



1 + c' 
et 



y = b- 



l + d' 



1 + c 



fu 
1 + d' 



(15) 



1 -r 2 1 r, 2 
a X A bPA 



e-^A-^PA^ (16) 



, (17) 

(18) 
(19) 
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Wd 



V(l + c)(l + d)' 
2 

V(i + c)(i + 4 e 



-t 2 /(l+c)-u 2 /(l+d) 



(20) 
(21) 



The photon subtracted states are just the difference 
between two squeezed states. For the total state, we can 
then gather terms to get 



W out (x A , PA ) =( /sq + ^isp 1 + (1 _ x) R sq 



1 (1 - x)-Rdis P + -R. 



dc 



(i-x)^ sq 



1 - Wd 
w 1 



R l-w 



R 



1 _l r 2 _ 12 
p a X A bPA 



1 



R l-w ttVo 7 ^ 



^ ^sq + ^disp 1 e - ^ (x A -r d ) 2 ~ £ (p.4 - «d ) 2 



R 



l-w d ttVo^ 7 



(22) 



PHYSICAL DESCRIPTION OF STATE BEFORE 
TRIGGER DETECTION 

In this section, the generalized variables of the co- 
variance matrix and displacement vector in the previ- 
ous section are related to experimental parameters. The 
two-mode Gaussian state before the trigger detection has 
temporal correlations given by the OPO output correla- 
tions ([5|-(|6|, by the trigger filtering and by losses. The 
APD click automatically defines a temporally localized 
single mode of the trigger beam, and based on the click 
time, we define a localized mode for the signal beam. Al- 
though this mode selection occurs after the actual photon 
subtraction, it is possible in the model to impose it on the 
initial Gaussian state, in order to get a time-independent 
covariance matrix and displacement vector. 



Covariance matrix 

The initial time-dependent covariance matrix is mod- 
ified by the tapping beam splitter, losses, and the mode 
selection. The handling of covariance matrix transforma- 
tions is simplified a little by using normal ordering - we 
do not have to take account of the vacuum contributions 
until the end. 

The initial two-mode time-dependent covariance ma- 
trix is 



•Topo^ — t f ) • = 



( :r OP o,ii: \ 

:r PO,2 2 : 



\ / 



(23) 



with 



:r po,ii: = 2(:Ax AAni (t)Ax AAni (t'):), 
:r po,22: = 2(:Ap A , ini (t)Ap AM (t'):). 



(24) 
(25) 



The tapping beam splitter, represented by the matrix 
V(T t ), transforms it to 



To obtain the state within the selected temporal mode, 
we integrate the quadrature variables over the relevant 
filter functions, /a f° r example, 



/oo 
fA(t)x A ,bs(t)dt. 
-oo 



(27) 



When calculating the elements of the covariance matrix 
r, we can move the integration outside the expectation 
values, e.g. 

:r 13 : = 2(:Ax A Ax B :} 

= 2(: Jf A (t)Ax Khs (t)dt Jf B (t')Ax B , hs (t')dt':) 

= [f A (t)f B (t') 2(:Ax Aths (t)Ax Bths (t'):)dtdt' 



f A (t)f B (t') :T bStl3 (t-t'):dtdt'. 



(28) 



Losses in both modes, as well as the trigger frequency 
filtering can be incorporated in the filter functions (see 
later). Therefore, the final matrix is 

:T:= JJf(t) T :T hs {t-t'):f{t')dtdt', (29) 

with f(t) = diag(f A (t),f A (t)J B (t),f B (t)), and, remem- 
bering to put back the vacuum: 



r = :r:+r vac . 



Displacement vector 



(30) 



With a displacement Db{P) = D B (|/3|e^ disp ) in mode 
B, the displacement vector is 

d = y/2\fi\ (0, 0, cos disp , sin (/> d i sp ) T . (31) 



The photon number in the mode due to this displacement 

Lisp 



is of course ndi sp = \fi\ 2 . On the other hand, from n = 



:r bs (t - t') := V(T t ) :T ini (t - t') : V(T t ) T . (26) 
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((: Ax 2 :) + (: Ap 2 :))/2, the photon number due to the 
squeezing can be obtained as 



1 



(=r: 



33 • 



:r 44 :). 



(32) 



The photon numbers are not directly experimentally 
accessible, but they are proportional to the detected APD 
click rates, so we can get \f3\ and the displacement vector 
in terms of these rates: 



^disp 



^sq 



' ^disp • r ; 



33 • 



44 • 



^sq 



R 



(33) 



sq 



d = 



' ^disp : : + : T44 : 



(0, 0, cos 0dis P , sin (j)disp) 7 



(34) 



Mode functions 



In the approach taken here, the filter functions play 
two different roles: 

1. To model the physical effect on the correlation ma- 
trix by losses and trigger filtering. 

2. To describe the actual observed temporal mode. 

#2 is the original meaning of the mode function, but #1 
can be included very nicely. 

For the signal (mode A), the observed mode is the 
temporal filter applied to the raw homodyne data in post- 
processing. This can be chosen freely, but we always use 
the experimentally most reasonable function which is 



Mi) 



-7l*l 



- K \t\ 



(35) 



with normalization Ma — 7 3 ft 3 (7+ft)/(7 4 +7 3 ft — 47 2 ^ 2 + 
7ft 3 + ft 4 ). This is the double-sided exponential from the 
OPO output correlations, smoothened by the frequency 
filtering of the trigger. Including the overall signal effi- 
ciency gives the filter function 



Ia{P) = y/VA^A(t). 



(36) 



The APD detection time is very short relative to the 
correlation time, so the temporal mode can be taken to 
be a delta function: 



The frequency filtering, however, has the effect that the 
photons arriving at the APD has been delayed relative 
to their correlated twins in the signal path. This can 
be taken into account by convoluting the mode function 
with the filter response, which gives (assuming the single 
Lorentzian filter) 

f B (t) = ^2^r] B e Kt , for t < and otherwise. (38) 

The overall trigger efficiency was also included in the 
filter function. 

Both filter functions assumes a photon click at time t = 
0. Due to the inclusion of losses, they are not normalized 
to 1. 



IDEAL SQUEEZED QUBIT 

When we compare the experimentally obtained or the- 
oretically predicted states with the ideal squeezed qubit, 
we calculate the overlap (fidelity) between the Wigner 
function of these states and the Wigner function of the 
ideal squeezed qubit 0. Given parameters r (squeezing 
parameter), and <j) (Bloch vector), the ideal squeezed 
qubit has the state ket 

|SQ> = cos 9 -S(r) |0> + sin °-S{r) |1> , (39) 

For the calculation of the Wigner function, W^Qix^p) = 
(27T)- 1 / e iyp (x - y/2 |SQ) (SQ| x + y/2)dy, we need the 
wave functions of the squeezed vacuum and the squeezed 
photon: 



(x\S(r)\0) = (ne») 



2r\-l/4-x 2 /2e 2 



(40) 



(x\S(r)\l) = (x\ 



1 



aS(r) |0> 

( d 
V^sinhr \ dx 



sinhr 
1 



(x\S(r)\0) 



-x 2 /2e 



3r 



(41) 



This gives the following expression for the squeezed 
qubit Wigner function: 



(37) 



W S q(x,p) 



cos t 



7T 



(i 



COS 



< 



,2r 



-2r 



) + V2( 



COS ( 



smc 



: ) sin O^j. 



(42) 
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